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N-mode coherence in collective neutrino oscillations 
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We study two-flavor neutrino oscillations in a homogeneous and isotropic ensemble under the 
influence of neutrino-neutrino interactions. For any density there exist forms of collective oscillations 
that show self-maintained coherence. They can be classified by a number A'' of linearly independent 
functions that describe all neutrino modes as linear superpositions. What is more, the dynamics is 
equivalent to another ensemble with the same effective density, consisting of A'' modes with discrete 
energies Ei with i = 1, . . . , A'^. We use this equivalence to derive the analytic solution for two-mode 
(bimodal) coherence, relevant for spectral-split formation in supernova neutrinos. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

In the early universe and in collapsing stars, the den- 
sity of neutrinos is so large that they produce significant 
refractive effects for each other. In a seminal paper 1], 
Pantaleone showed that this is not just a correction to 
the usual matter effect 0-111 , but causes qualitatively 
new phenomena. They originate from flavor off-diagonal 
refraction caused by neutrino oscillations which there- 
fore feed back on themselves. One consequence is self- 
induced flavor conversion even for a very small mixing 
angle, caused by an instability of the interacting ensem- 
ble. When the density decreases slowly, for example in 
a supernova as a function of radius, this instability leads 
to flavor swaps between neutrino spectra in sharp en- 
ergy intervals, causing conspicuous "spectral splits." The 
body of literature on collective oscillations has imme nsely 
grown, so we cite only some theoretical key papers |5l-l26| , 
whereas we refer to a recent review (27j for the torrent of 
activities in the area of supernova neutrinos. 

One key feature of collective oscillations is that, given 
suitable initial conditions, the entire ensemble evolves 
in a highly correlated way (self- maintained coherence). 
We make this concept more precise, if only for the sim- 
plest case of a homogeneous and isotropic gas with two 
mixed flavors and ignoring ordinary matter. Moreover, 
we mostly consider fixed density, although we use moti- 
vations related to adiabatic changes of density. 

In the absence of neutrino-neutrino interactions, every 
mode evolves independently of the others. The ensemble 
may be described by two-spinors in flavor space V'-E(i), 
occupation- number matrices QE{t), or the corresponding 
polarization vectors E{t) — effectively there are three 
important numbers for every mode: the amplitudes of 
the two flavor components and their relative phase. We 
always use polarization vectors and label the modes with 
their vacuum oscillation frequency to = /S.m? where 
negative uj denote antineutrinos. Vacuum oscillations are 
described by the polarization vector of a given mode pre- 
cessing with frequency w around a direction given by a 
unit vector B, the mass direction in flavor space. If neu- 
trinos begin in weak-interaction eigenstates, all P;^ are 
initially tilted relative to B by twice the vacuum mix- 



ing angle (Fig. [T]). Each precesses with a different 
frequency and so eventually the modes distribute them- 
selves uniformly on the precession cone centered on B, 
approaching complete kinematical decoherence. 

Neutrino-neutrino interactions provide an additional 
force in terms of the global P = J duj P^ , leading to the 
equations of motion (EoMs) 



P^ = {loB + ^P) X P, 



(1) 



Here fj, ~ v^Gpn^ is a typical interaction energy, its 
exact value depending on the normalization of the polar- 
ization vectors. In the limit /i ^ oo, and if initially all P^^ 
are aligned in the flavor direction, they "stick together" 
and precess around B with the common frequency (w). 
These "synchronized oscillations" are the simplest case 
of self- maintained coherence |Tl|, [l^] ■ 

Starting with such a configuration, can slowly de- 
crease as in the expanding universe or with distance from 
a supernova. The ensemble continues to oscillate coher- 
ently, but the Poj slowly separate and fan out in a plane 
that precesses around B. Such a pure precession mode 
can exist for any strength of /j, [l^ . If ^ decreases adi- 
abatically to zero, the final configuration is that all P^^ 
with UJ above some frequency Wgpiit point in the positive 




FIG. 1: Precession of a polarization vector P^^ with frequency 
(jj around the z-direction (the mass direction in flavor space). 
The initial orientation is in the weak-interaction direction, 
tilted by twice the mixing angle relative to the ^-direction. 
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B direction, the others in the opposite direction, forming 
a "spectral split." 

More complicated forms of coherent motion occur in 
the context of multiple split formation [23]. A generic 
example is a gas of neutrinos and antineutrinos of a given 
flavor with a thermal distribution. Figure [2] shows a 
Fermi-Dirac distribution as a function of the re-scaled 
frequency uj = T / E. In this way u; is a dimensionless 
variable of order unity. The neutrino flux streaming from 
a supernova core contains an excess of Vf, over Pe because 
the trapped electron-lepton number needs to be carried 
away. We mimic this situation schematically by assum- 
ing a small degeneracy parameter 77 = 0.2, corresponding 
to ~ 0.70 The region around lu — corresponds 
to the high-i? tail and the high-w tails (small E) fall off 
as w"^. The mixing angle is taken to be very small and 
so the mass and flavor directions are almost identical. 
Notice that P„ pointing up means one flavor, pointing 
down the other. In the flavor-isospin convention [16, 23], 
the interpretation is reversed for antineutrinos, explain- 
ing that the spectrum in Fig. [5] is negative for w < 0. 




-0.5 0.0 0.5 
Frequency w 



FIG. 2: Neutrino spectrum (z-component of P^j) in the vari- 
able UJ = T / E, assuming that initially only one flavor is popu- 
lated with a Fermi-Dirac distribution (degeneracy parameter 
77 = 0.2). The integral over positive uj is normalized to unity. 



Without v-v interactions only small-amplitude vacuum 
oscillations occur, whereas for non-zero fi the system is 
unstable. The evolution is nicely illustrated by show- 
ing the instantaneous "swap factor" Zi^{t), i.e. the w- 
dependent factor by which the original spectrum must 
be multiplied to obtain the instantaneous z-component 
of P<,(0, i.e. P^{t) = z^{t)P^{Q) dl . The entire spec- 
trum oscillates coherently. In the top panel of Fig. [3] we 
show numerical results — 10) for the maximum and 
minimum z^ during an oscillation period as well as their 
average. 

In this example and in any unstable spectrum, z^^ (t) is 
initially a Lorentzian of a certain width A that oscillates 
like an inverted gravitational pendulum with natural fre- 
quency equal to the same A [23| . When ji decreases adia- 
batically, z^ changes, the oscillation amplitude decreases, 
and at = takes on square- well shape: a spectral swap 
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FIG. 3: Maximum, minimum and average swap factor for 
the thermal spectrum of Fig. [2] if the interaction strength 
begins with fj, = IQ and then decreases adiabatically to 0. 
Snapshots for fj, = 10, 3, 1, 0.3 and (top to bottom). The 
initial spectrum is overlaid using a vertically compressed scale 
relative to Fig. [S] 
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with two sharp sphts at its edges has formed. Several 
snapshots are shown in Fig. [31 

These examples show "coherence" in some intuitive 
sense, but what precisely does this mean? We here define 
coherent motion as a form of evolution where neighboring 
modes do not separate from each other. Kinematical de- 
coherence means that different Puj{t) within a frequency 
interval Slo eventually separate by a large amount, no 
matter how small Suj. If we coarse-grain the ensemble 
with a bin width 600, the coarse-grained (Puj) becomes 
shorter than as time goes on [2^. On the other 
hand, coherence means that the coarse-grained ensemble 
remains identical with the original one, no matter how 
much time has passed. Pure precession is a trivial exam- 
ple because the polarization vectors remain static relative 
to each other, whereas in the example of Fig. [31 distant 
Puj{t) move wildly relative to each other, yet neighboring 
ones do not separate as time goes on. 

Independent evolution of neighboring modes, no mat- 
ter how close their frequencies, implies that every Puj{t) 
is linearly independent. For vacuum oscillations this is 
obvious because each Puj{t) is a harmonic function of 
frequency uj. On the other hand, kinematical decoher- 
ence can not occur if every P^it) is a linear combination 
of a small number N of independent functions. In this 
case the system performs what we call A'^-mode coherent 
oscillations. The infinitely dimensional function space re- 
quired for kinematical decoherence has reduced to a small 
number of dimensions. 

What is more, the same function space is spanned by 
the solutions 3i{t) of another system {i ^ 1, . . . ,N) with 
vacuum frequencies uji and the same effective density. So 
we consider the EoMs of a new system 

j, = {uj,B + n3) X J, . (2) 

The functions 3i{t) can be found such that 

m - j(o , (3) 

meaning that both systems are dynamically equivalent. 
The original modes Pu{t) and the linearly independent 
"carrier modes" 3i{t) span the same function space: each 
Pui{t) is a certain linear combination of the Ji(i). Notice 
that usually the carrier modes can not be interpreted as 
a coarse-grained representation of the original ensemble, 
but rather are more abstract. 

The idea that coherence in collective neutrino oscilla- 
tions corresponds to linear dependence among different 
modes and that the evolution is dynamically equivalent 
to a small ensemble of carrier modes provides a simple 
picture for self-maintained coherence. Our primary inter- 
est is to use this approach for the construction of explicit 
classes of solutions, notably for bimodal coherence. 

We begin our study in Sec. [ill with the equations of 
motion, the idea of A^-mode coherence, and the concept 
of carrier modes. Then we proceed to the explicit cases 
of single-mode, two-mode and multi-mode coherence in 
Sees. [mHYl and conclude in Sec. [VD 



II. N-MODE COHERENCE 
A. Equations of motion 

The two-fiavor evolution of a homogeneous and isotro- 
pic neutrino gas can be described by flavor polarization 
vectors Ptj(i), where lu = Am^ /2E is the vacuum oscilla- 
tion frequency. P^^ pointing up means one flavor, point- 
ing down the other. In the flavor-isospin convention jT6| . 
the interpretation is reversed for antineutrinos (uj < 0). 
The flavor evolution under the influence of vacuum os- 
cillations and neutrino-neutrino interactions is given by 
Eq. ([T|. This precession equation has been stated many 
times[i, 11, i, fill [H E, ii, [13, most recently in great 
detail in Ref. [26j, and we use it here without further ado. 

A few points still deserve explicit mention. The polar- 
ization vectors derive from an ensemble average and thus 
are classical quantities. Our EoMs rely on a mean-field 
assumption: a given neutrino is supposed to feel only 
the average effect of all others, ignoring fiuctuations or 
quantum correlations. (The quantum to classical transi- 
tion was studied in Refs. [H [11.) The length of P^(i) 
is not fixed by a quantization condition, but chosen at 
convenience. The value of fi depends on this choice. 

The EoMs conserve energy and derive from the classi- 
cal Hamiltonian 

H = J LuB-P^dLo+^-^P\ (4) 

Here, H and all P^ are functions of the canonical coordi- 
nates and momenta. The evolution of any function F on 
phase space is given by the Poisson bracket F = [F, H] . 
It is crucial that the Pi^ play the role of classical spins, 
or rather isospins. They obey the angular- momentum 
Poisson bracket [Puj,x, P^' ,y] = S{uj — w') Pi^^z and cychc 
permutations. Notice that [ • , • ] is not a commutator — 
classical variables commute in the sense PxPy — PyPx- 
The overall P is the total angular momentum (in flavor 
space). Collective oscillations in the mean- field approxi- 
mation are thus equivalent to a set of interacting classical 
spins. However, the simplicity of our Hamiltonian is de- 
ceiving because it encapsulates all the complications of 
collective fiavor oscillations. 



B. Rotating frames 

Another conserved quantity is P projected on B, 
equivalent to angular-momentum conservation in the 
symmetry direction. Therefore, we always work in the 
mass basis and use coordinates where the z-direction co- 
incides with B (Fig.[T]). Usually all P^j begin in a weak- 
interaction eigenstate and thus in a common direction 
that is tilted relative to B by twice the mixing angle. 

The conservation of Pz = B-P implies that the internal 
motion of the ensemble, up to overall precession, does not 
depend on P^. Depending on convenience, we can study 
our system from the perspective of a frame co-rotating 
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with some frequency [l5| . We may choose Wc = /iPz and 
without loss of generaUty study the equivalent EoMs 



(5) 



where Pj_ is the part of P transverse to B. We will fre- 
quently encounter the degeneracy between a shift of the 
neutrino w spectrum by some amount Wc and ^iPz- On 
the Hamiltonian level, such shifts amount to subtracting 
the conserved quantity WcB • P. 

Shifting a spectrum such as Fig. [2] by an amount ujc 
modifies the interpretation because negative-frequency 
modes (anti-neutrinos) become positive-frequency modes 
or vice versa. A mode that used to be occupied by Df, 
is then interpreted as being occupied, for example, by 
v^. However, this modified interpretation does not affect 
the abstract dynamics, and we can always shift back at 
the very end. The possibility to shift anti-neutrino and 
neutrino modes seamlessly into each other is the main 
advantage of the fiavor isospin convention. 



C. Linear dependence of different modes 

Collective neutrino oscillations consist of strong corre- 
lations among all modes and we have advanced that this 
means that the F^j (t) depend linearly on a small number 
N of independent functions. To be more precise, we use 
N to denote the number of nontrivial functions because 
in addition the constant B will be seen to appear, so 
Pui{t) actually depends on + 1 functions. 

It is illuminating to diagnose N in numerical examples. 
The ensemble is represented by a large discrete set P^ (t) 
with i = I, . . . ,n and n 3> 1. To find the number of lin- 
early independent functions on a time interval ti < t < t2 
we calculate the Gram matrix Gij = J*^ dt P^ {t) ■ Vj {t) . 
Its rank reveals the number of linearly independent func- 
tions. In reality, one finds N -\-\ large eigenvalues of Gij 
and the rest very much smaller. They can be made ever 
smaller by increasing the degree of adiabaticity, i.e. slow- 
ing down ^,{t) and decreasing the mixing angle, i.e. mak- 
ing the initial polarization vectors more collinear with 
B. Within constraints of numerical accuracy, the chosen 
time interval is irrelevant to find -I- 1 if the P^ [t] are 
analytic functions, but of course infiuences the numerical 
realization of the Gram matrix. 

We have always found the expected A^. The forma- 
tion of a single spectral swap as in Fig. [3] corresponds 
to bimodal coherence at any stage, i.e. the Gram matrix 
was found to have rank 3. The evolution from a double- 
crossed spectrum towards two swaps [23 | corresponds to 
four- mode coherence (rank 5), and so forth. 



D. Carrier modes 

If every ^^{1) is a linear combination of N inde- 
pendent functions and the constant B we may choose 



any set V^-it) with i = 1,...,N as a basis if they 
are not accidentally degenerate. We then have P{t) — 

J dujP^{t) ~ cqB + X^fci aiPwi(i) as a unique lin- 
ear combination. We define Ji{t) = UiPujiit), obeying 

Ji = {uji'B + ao/iB -I- fj,3) X where J = X^i^i Jj- & 
rotating frame we can transform away ag/iB, or rather, 
we may choose another set Pi^'{t) where u'q —0. 

So we can always find a set of linearly independent 
functions 3i{t) with i = 1, . . . ,N that obey Eq. ([2|), ful- 
fill the matching condition Eq. ([3]), and together with 
B span the same function space as the original ensem- 
ble Pcj(i). The internal dynamics of our system is the 
same if we study the co-rotating EoMs of Eq. ([5]). In this 
case the matching condition is Pj.(^) — J±{t) and the 
complication about the oq coefficient disappears. 

We can express each P,^ (t) of the original ensemble as 
a linear combination of the chosen "carrier modes" Ji (t) 



N 



(6) 



i=l 



This function must obey the original EoM of Eq. ([ij. 
Inserting Eq. ^ on the l.h.s. of Eq. (P) and using Eq. ^ 
yields J2i=i ^i('^iB -I- /xJ) x Ji. On the r.h.s. of Eq. ^ 
we use P = J and Eq. ^ such that —bofiB x /iJ + (a;B + 

^J) X J2iLi ^jJj- With J = X]t=i "^i fi'^st term and 

after collecting everything we find 



N 



= B X ^[bi (uj - uji) - ban\Jiit) . 



(7) 



The functions 3i{t) are linearly independent by construc- 
tion and not proportional to B. Therefore, we find the 
coefficient bi — b^^/ {lj — uJi) and thus 



Pu{t)=ba 



N 



■J.(t) 



(8) 



The required linear combination is unique up to an over- 
all proportionality factor. (This linear combination is 
an example of a more general transformation brieffy dis- 
cussed in Appendix [Cl) 

The expression, Eq. ([81), is singular at the carrier fre- 
quencies. However, it is only the direction that matters, 
so we introduce the transformed carrier spectrum 



N 



1=1 



N 



1=1 



where the overbar has nothing to do with antiparticles. 
Introducing the spectrum g^) with \Pui\ = \guj\, the origi- 
nal ensemble in terms of the carrier modes is 



P - o — 



(10) 



where 
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We usually consider initial conditions such that the P^i 
are almost coUinear with B and we ask for the probabil- 
ity for a neutrino mode lu to stay in its original flavor, 
perhaps after an adiabatic change of This information 
is contained in the z-component of J^j and we call 



'Jlj 



(11) 



the time-dependent "swap factor." 

Our discussion is motivated by the spectral swaps and 
splits that form when the effective density encoded in 
/i slowly decreases to zero. Since we are dealing with a 
Hamiltonian system, we expect the adiabatic change of 
the parameter fj, to deform the solution, but that it re- 
mains coherent if it was initially coherent. For /x — ;> 
we have J;^ oc B, implying that in the end all polariza- 
tion vectors are collinear with B. Conceivably one could 
define A'^-mode coherence by the very property of devel- 
oping N spectral splits when fi decreases adiabatically to 
zero. 

Our main interest is to use these ideas to construct 
explicit classes of A^-mode coherent solutions, notably 
the most general bimodal solution. For a set of functions 
3i{t) fulfilling Eq. the matching condition Eq. ([3]) is 



3{t) 



dujgu. 



(12) 



For a given spectrum and interaction strength /i this 
condition restricts possible sets of carrier modes to pro- 
vide an TV-mode coherent solution. 



III. SINGLE-MODE COHERENCE 

A first trivial case is what we may call zero-mode co- 
herence, where all P^j are exactly aligned with B and 
therefore stay that way. This configuration also appears 
as a limiting case of other forms of oscillation. 

The first non-trivial case is single-mode coherence, 
identical with pure precession 113 ■ The single car- 
rier mode 3{t) precesses around B with some co- rotation 
frequency uJc and obeys J = x J where the nonlinear 
term /it J x J = has dropped out. According to Eq. (|9]) 
the transformed spectrum is J^j = (w — Wc) B -|- /iJ. The 
matching condition P = J then implies 



{lo - Wc) B + ^P 



(13) 



and thus = (w — Wc -I- nPzY' + {pP±)^, where Pj^ 
denotes the part of P transverse to B. The matching 
condition Eq. (|12p is, in agreement with Ref. (22| . 



Pz ^ J dujg^ 
Pi = / duQ„, 



UJ -OJC+ nPz 



^{uj-UJ,+llPzf + {lJiPi_Y ' 
1^ 

^{u-u, + iiPzY + {iiPA_Y ■ 



(14) 
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FIG. 4: Swap factor for pure precession with k, = Q, 0.5 and 1 
according to Eq. (|16p . 



We have used the fact that all P;^ lie in the plane spanned 
by B and P. 

We could have studied our system from the perspec- 
tive of a rotating frame as discussed in Sec. IIIBI In the 
above expressions we recognize that nothing changes if 
we absorb ^.Pz in the definition of uic- 

For a given spectrum g^ we can use these equations 
in different ways. If Pz is fixed by an initial condition 
we may ask for the value of P±_ and lOc corresponding 
to a specified value of /i. In particular, for /i — >■ we 
can find the split frequency for a given Pz- For explicit 
solutions in this sense and further discussions we refer to 
the literature [l7l[22|. 

We may also use these conditions to construct all pos- 
sible pure precession solutions for a given g^^ . To this end 
we introduce the parameters k = /LtPj_ and 7 = — ^Pz , 
leading to the general pure precession solution in the form 
of a two-parameter family depending on 7 and k, 



Plo , z — gu. 



Puj,i. = gu. 



w — 7 



v/(a; - 7)2 + ^2 

K 

V(W - 7)2 -I- ^2 



(15) 



In addition one needs to specify the precession phase and 
thus the instantaneous orientation of the comoving plane. 
Integrating over du yields Pz and and the correspond- 
ing /J, = k/ P_l and Wc = 7 + ^jiPz ■ 

In terms of the parameters k, and 7, the swap factor 
for pure precession is 



w — 7 



\J{lO — 7)2 + k2 



(16) 



Some examples are shown in Fig. |4l The step occurs at 
frequency 7 and the transition region has width k. 

For K — > the swap factor is identical with the sign 
function z^^ — sgn(aj — 7). In this case the first relation 
in Eq. (fT4)) has the form Pz = J dui g^ sgn(a; — 7) and can 
be used to find 7 for a given P^. The second equation 
is 1 = J duj guj+j id/^/uj"^ + k2 and can be used to find k 
explicitly in the limit where fj, is very small. In this limit, 
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the integral is peaked around a; = and we can approx- 
imately pull out g.y. The remaining integral diverges for 
ljj ±cx), so we cut it off at some frequencies ±a, leading 
to 1 ^ 2/i(7^ log(2a/K) and implying that for small 



,-l/(2MffT) 



(17) 



In other words, k becomes exponentially small for small 
/i or small g^. 



IV. BIMODAL COHERENCE 

A. Matching conditions 

Turning to bimodal coherence, we consider two carrier 
modes Ji,2(i) with frequencies u;i,2 fulfilling the preces- 
sion equation, Eq. According to Eq. ^ the trans- 
formed carrier spectrum is 

(a; - ClIi)(w — a;2) B 

+ /i (w — a;2) Ji + ^ (cj — wi) J2 . (18) 

Two interacting polarization vectors are dynamically 
equivalent to a gyroscopic pendulum (Appendix |X| . To 
use this equivalence we write wi = Wc— /3 and W2 = ^c+P, 
where Wc is a suitable co-rotation frequency such that the 
two carrier modes have equal but opposite frequencies 
±/3. To simplify notation we implement loc by shifting 
the spectrum, i.e. instead of g^, we use g^+Lo^ in the inte- 
gral for the matching conditions. 

From Appendix VK\ we borrow Q = Ji — J2 + (/3//-i) B 
for the pendulum's radius vector and find 



Jt^ =w^B + /ia;J + /i/?Q, 



(19) 



where J = Ji -I- J2 is the total angular momentum. To 
spell out the matching conditions we need a unit vector 
in the J^j direction and thus need 



- 2/iw^B • J + 2i?l3uj Q • J . 



(20) 



From Eq. (jA6l) we take the natural frequency = /3/i|Q|, 
so the second term is A"*. The bracket in the third term 
is the energy E of the pendulum. In the first term on the 
second line we have B • J = , a conserved quantity, and 
in the final term Q • J = |Q| 5 with S the spin (conserved 
angular momentum projection on the pendulum axis). 
Ordering by powers of lj we find 



3l^uj^ + 2n {io^J, + Lu^E + ioX'^S) + A-* 



(21) 



As expected, this quantity is conserved, i.e. it does not 
depend on the oscillation phase. 

We also need the individual components of J^j. We 
select directions in the plane transverse to B that move 
with the system and we call the direction along Jj_ the co- 
moving ?/"direction, where is the part of J transverse 



to B. The comoving a;-direction is then orthogonal and 
hence collinear with B x J. We project out the three com- 
ponents by taking the scalar product of J^j with B x J, 
J 1 = J — JrB and B and find 



Jui.T. 



J. 



(B X J) 



Ju.z = oj'^ + ^ii^Jz -t- Q • B . 



(22) 




On the l.h.s. of the matching condition Eq. (IT^ we have 
in these comoving coordinates J = (0, Jj^, Jz) so that 



(23) 



In the first line we have divided out Ju>,x because it does 
not depend on w, and using this condition then simplifies 
the second line, and both are used to find the third. So we 
could find the matching conditions without even spelling 
out the scalar products in Eq. ([2^ . 

Every solution is given in terms of the four pendulum 
parameters E, S, and A, although we have not spelled 
out Eq. ([22I) in terms of these conserved quantities. In 
addition we have the shift frequency Wc as a fifth param- 
eter. However, we will see that ujc and jiJz appear only 
combined, similar to the pure precession case, so actu- 
ally there are only four independent parameters. The 
dynamics is given in terms of cos 5, the angle between 
B and Q, that follows the pendulum differential equa- 
tion of Eq. (|B5p . In addition we can calculate the pre- 
cession motion following the usual pendulum treatment 
(Appendix [B|. 

Given the pendulum parameters and the EoM for cos 9 
we have all the information needed to describe the mo- 
tion. However, these parameters do not allow us to deter- 
mine the original carrier modes in a unique way — there 
are many solutions. Notice, in particular, that the fre- 
quency P and the length of Q are degenerate parameters. 
Two interacting polarization vectors are equivalent to a 
pendulum in flavor space, but the reverse mapping is 
not unique. In particular, the carrier frequencies are not 
unique. In certain limits, however, a particular choice 
may be most convenient. 



B. Swap parameters 

While the pendulum parameters are well matched to 
the underlying dynamics, they are not very intuitive 
when studying cases such as the one shown in Fig. [3) To 
find a more suitable parametrization we first spell out 
the swap factor Zi^{t) = Juj,z/ Ju> explicitly. 



Zu:{t) = 



LtP' + yucj + A^ COS 



^a;4 -f 2/^ (a;3 Jz -f uj'^E ^ ujX^S) -f A"* 



(24) 
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Inspecting the numerical results of Fig. [3l the analytic 
swap factor for pure precession of Eq. motivates the 
following ansatz for the average swap factor 



(w - 7i) 



- 72) 



(25) 



where 71^2 give the location of the two steps and ki_2 
their respective widths. 

Comparing the coefficients of the lo polynomial in the 
denominators of Eqs. (1^^ and provides 



A* 



ill + ^^Diii + 4) , 
-(71 + 72) , 

(7i+72)^ + 27i72 + K? + Ki 



71(7 



■72(7? 



A2 



(26) 



A further simplification is achieved if we express the time 
variation of the swap factor in terms of a new dimension- 
less variable 



cos 9{t) - 7172 



K1K2 



implying 



(g; - 71) (^ - 72) + KiK2u{t) 



(27) 



(28) 



Inserting u{t) in the EoM for cos6{t) of Eq. (jBSp and 
using the parameter mapping of Eq. (|26p provides 



= (1 - 1*2) [(^^ _ ^2)2 + (ki - K2f + 2kiK2(1 " «)] , 

(29) 

where the range of motion is— 1<m<-|-1. Therefore, 
the average is w = and the average swap factor is indeed 
given by the proposed expression Eq. ([25l) . 

The structure of Eq. ([29]) suggests the interpretation 
u — cos?? with d{t) some abstract angle variable. In this 
way \ — = SIT? § and one finds 

= (71 - 72)^ + (ki - ^^2? + 2niK2(l - cos I?) . (30) 

This angle always advances in the same direction with a 
modulated velocity. The EoM is that of a gravitational 
pendulum with a speed so large that it moves through 
the inverted position. 

We show examples for the maximal, minimal and av- 
erage swap factors, corresponding to u = -1-1, and —1 
in Fig. [5l where 71 = —1, 72 = and ki — 0.1 are held 
fixed and we show different cases for H2 — 4, 2, 1, 0.5 and 
0.1 (top to bottom). 

We realize that in Eqs. ([25)1 and (1^^ the parameters 71 
and 72 always appear either as 71 — 72 or as w — 71^2- So 
we can go to a frame rotating with some frequency Uc by 
shifting u and 71^2 by this amount and the solution will 
be the same. Since —(71 +72)//^ has the interpretation of 




-2 2 

Frequency oj 



FIG. 5: Maximum, minimum and average swap factor accord- 
ing to Eq. (f28)) using 71 = —1, 72 — -1-1, ki = 0.1, and K2 ~ 4, 
2, 1, 0.5 and 0.1 (top to bottom). 
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Jz this means that we can trade for a co-rotation fre- 
quency Wc • Once more this is the same effect discussed in 
Sec. Ill B[ Instead of fixing = it is enough to match 
P_L = Jx and trade the mismatch of z-components for 
a suitable co-rotation frequency Wc- Here, in particular, 
we can adjust Jz such that uJc = 0. 

Spelling out the matching conditions in terms of the 
new parameters, we find 





P. 



= du! ■ 



1 



- 7i)(w - 72) 



(31) 



o 
CO 



CO 




where 



(32) 



FIG. 6: Swap factor for pure precession as derived from the 
gyroscopic pendulum picture. The chosen parameters are 
71 — —1, 72 = -1-1, Ki = and K2 = 0.5. 



For a given spectrum, the first equation establishes a con- 
straint among the four parameters ki^2 and 71,2- The 
second allows us to calculate the required /i for this solu- 
tion, and the third provides the required Pz ■ We can also 
use the first and second condition to simplify the third. 

It is also of interest to spell out explicitly all three 
components of J;^ of Eq. ([2^. For Jj^ we find 



and 



(33) 



up to an overall precession. The swap factor, no longer 
depending on time, is 



sgn(a; - 71) 



(w - 72) 



\/(w - 72)2 + k: 



(35) 



where we have used ki = 0. In other words, we find 
the pure-precession swap factor of Eq. (jl6l) . augmented 
with a step function centered on a frequency 71 (Fig. [6]). 
Therefore, our previous result is equal to the present one 
for the choice 71 — >■ —00. 



Ju),y — 



\/{l2 - 71)^ + «S 



K1K2VI 



{2uj - 71 - 72) kI + (72 - 1i){kI - kI) 
{oj - 7i)(cj - 72) + ■ 



(34) 



Divide these components by J^o of Eq. ([32]) to find the 
unit vector ^ujI Juj- 

For — 00 < w < +00 this vector traces a closed curve 
on the unit sphere, where the points w = ±00 are at the 
north pole. This closed curve moves as a function of i9(i). 
Notice that for u — ±1, corresponding to the highest 
and lowest swap-factor curves in Fig. [31 the comoving 
x-component vanishes for all w, implying that all lie 
in the plane spanned by B and P. The case u — —\ 
corresponds to the lower swap factor, tracing out all z- 
values between —1 and -1-1, i.e. the unit vectors trace out 
the great circle spanned by B and P. 



C. Pure precession limit: k;iK2 = 

Bimodal solutions are described by two carrier modes, 
which in turn are equivalent to a gyroscopic pendulum. 
One possible form of motion is pure precession: single- 
mode coherence appears as a limiting form of two-mode 
coherence. In terms of our swap parameters, this case is 
described by ki = or K2 = where the solution is static 



D. Pure nutation limit: k\ = K2 

Another limiting case is when the swap factor is sym- 
metric around uJc — ^(71 + 72), i-e. when ki = K2 = k. 
In a frame co-rotating with Wc this solution is equivalent 
to a plane pendulum where Jz = S = and was found 
previously [23]. Shifting uj and 71,2 by ojc reproduces the 
pure nutation solution 



Vw* + 2a;^A^Cm + A" 



tjA^2(cm - c) 

^2 + A^c 

' (36) 

where s = sin 9, c = cos0, and Cm = cos^min- In contrast 
with Ref. [2^ we here count the zenith angle 9 from the 
north pole. In the co-rotating frame we have 71 2 = ±7 
and therefore the natural pendulum frequency is A" — 

+ 7^, whereas the plane pendulum's highest point is 
given by Cm = (k^ - 7^)/('«^ + 7^)- 

If the system has been set up with all polarization vec- 
tors initially almost aligned with B, the initial swap fac- 
tor is unity. The initial swap parameters are 71 = 72 = 
(jjc and Ki = K2 = K and are solutions of [23[ 




1 



1 



UJ — UJc 
{uj — UJc)^ + 



(37) 
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For the example of Fig. |3] with the initial value /i = 10 
one finds — —1.604 and k = 2.179. 



E. Two pure precessions: ki,2 <C I72 — 71 1 

When we begin with a bimodal system and the den- 
sity decreases adiabatically to zero, the end state shows a 
spectral swap with two splits at its edges (Fig. [3]). When 
II has become sufficiently small, the two quasi-step func- 
tions are well separated and the overall swap factor looks 
like the product of two pure precessions. This observation 
was the very motivation for writing the general bimodal 
solution in terms of the swap parameters 71.2 and K1.2 
instead of the pendulum parameters A, i?, Jz and S. 

In the present limit of ki,2 ^ I72 — 7i|, the two split 
regions are essentially independent. The carrier modes 
Ji,2(i) precess essentially freely with frequencies 71,2 ■ 
To study this case it is easiest to represent the x-y- 
components of all polarization vectors as complex num- 
bers of the form + i , so the transverse components 
of the carrier modes can be chosen as Kie''''^* and ^26''''^* 
with K\ — /ijJi^j^l and K2 = m|J2.±|- The abstract angle 
characterizing the overall motion is here 1? = ±(71 — 72 )t 
because in the present limit — (71 — 72)^. 

In the lower panels of Fig. |3]we see that the swap factor 
is asymmetric and steeper in the region where the step 
falls into a spectral region where is smaller. This is 
explained by Eq. ([T7|) where we have shown that for a 
pure precession with very small the width of the step 
is exponentially smaller for smaller g,^ in the step region. 

As decreases adiabatically in a case like Fig. |31 the 
system begins in a state of pure nutation and ends in a 
state of two independent pure precessions. These limiting 
forms and all intermediate cases are described by our 
four-parameter analytic solution. 



F. Numerical examples 



0) 1 







to 

1-1 



-2 



1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


















~ . . . 1 
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FIG. 7: Numerically determined swap parameters 71,2 and 
Ki,2 for the example of Fig. [3] as a function of /i that adiabat- 
ically decreases from 10 to 0. 




FIG. 8: Same as Fig. [T] now with /i being decreased from 
/io = 10 to zero as before, or increased from /xq = 10 to 
infinity. 



function of in Fig. |8l assuming the initial value is /iq = 
10. The solutions for < ^ < fiQ are from Fig. [T] We 
identify the parameters such that 71 < 72. 



To test if our analytic solution indeed corresponds to 
numerical examples such as Fig. [31 we extract the pa- 
rameters 7i_2 and ki,2 as fi slowly decreases (Fig. [71). 
Comparing the analytic swap factor with the numerical 
one yields perfect agreement. 

Usually we begin with P° almost aligned with B and a 
chosen /zq. Then /i decreases from fj,o to 0, going through 
different solutions as in Fig. [31 where remains con- 
served. If in this example we were to begin with /io = 3 
instead of 10, the /i = 3 solution would be the one with a 
Lorentzian pattern. So for the same gi^, Pz and fi there 
exist different solutions that we can construct numeri- 
cally by adiabatic deformations. 

We may also begin with a certain /j,o and then increase 
the density as during supernova collapse. The initial pure 
nutation is then deformed to another bimodal solution 
within our general four-parameter family. For our usual 
example of Fig. [31 we show the swap parameters as a 



V. MULTI-MODE COHERENCE 

A. Four-mode coherence 

More complicated forms of motion occur if the spec- 
trum has two or more independent instabilities, typically 
for multi-crossed spectra [23. Such spectra arise natu- 
rally if we augment the Fermi-Dirac spectrum of a single 
species (Fig. [^J with an initial population of the other 
species, but with a larger temperature and smaller flux. 
This mimics neutrinos streaming from a supernova core 
where a larger flux of Ve than De stream away and the 
other flavors and Vt have smaller fluxes, larger av- 
erage energies, and no asymmetry. So if we add to the 
asymmetric spectrum of Fig. [21a symmetric component of 
the other flavor with a 1.25 times larger T and 0.8 times 
smaller number density, the difference spectrum relevant 
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-1.0 -0.5 0.0 0.5 1.0 
Frequency co = T/E 



FIG. 9: Difference spectrum between two flavors. One is a 
Fermi-Dirac distribution (77 = 0.2) as in Fig. (2] the other 
a nondegenerate spectrum with 1.25 larger T and 0.8 times 
smaller f density. The integral over positive uj of the degen- 
erate flavor is normalized to unity. Overlaid is the flnal swap 
factor after /x has decreased adiabatically from 3 to 0. 



for flavor oscillations is shown in Fig. [9l 

This spectrum has two zero crossings with positive 
slope, allowing for two independent instabilities, i.e. 
Eq. ([37]) can have two solutions [23| . If /x decreases adia- 
batically from some value fj,o to 0, one finds two spectral 
swaps and four concomitant splits. As an example we 
show the final swap factor for /xq = 3 in Fig. [S] 

This case is an example for four-mode coherence. The 
complicated motion of all P^^ is equivalent to four carrier 
modes as confirmed by numerical tests using the Gram 
matrix. When fi has become very small and the swaps are 
almost complete, we have the now-familiar pattern of four 
independent pure precessions in the four split regions. In 
the final small-/i limit, the swap factor is 

z^^T] , '^"^^ (38) 

corresponding to four freely precessing carrier modes. 

Equation (|9]) implies that is now a polynomial with 
leading term w^, so the above representation for z^^ is 
general, except that in the numerator complicated time- 
dependent terms appear when we are not in the small-/x 
limit. One can group the four carrier modes in pairs and 
understand the overall motion as two gyroscopic pendu- 
lums interacting by a dipole force. Apart from an overall 
precession, the time dependence is then described by two 
nutation angles and one relative precession angle. 

If we choose some value for /xq, it is not assured that 
the system indeed has two instabilities — it can have only 
one or none. Therefore, as ^ adiabatically decreases, the 
system can at first show bimodal coherence and at some 
critical /i- value the second unstable mode kicks in, taking 
the system to a state of four-mode coherence. 



B. Stability issues 

It is generic that a system of lower modality can de- 
velop higher modality by an instability. Usually we be- 
gin with all polarization vectors almost aligned with B. 
If they were perfectly aligned, we would have zero-mode 
coherence for any spectrum g^^, but the smallest distur- 
bance allows the unstable modes to grow exponentially. 
We implement this disturbance in the form of a small 
mixing angle, i.e. a small mismatch between the initial 
condition and exact zero-mode coherence. The solutions 
(wc, k) of Eq. (P7|) identify the unstable modes, each one 
corresponding to a contribution of 2 to the total modality 
N relevant for the given g^^ and /i. 

We can also prepare a state of pure precession. With 
an arbitrary spectrum g^ and parameters 7 and k, the 
initial condition is defined by Eq. ([T5|) . allowing us to 
calculate the required /x. Henceforth the system evolves 
in single-mode coherent motion unless it is unstable. De- 
pending on g^ and with the slightest disturbance it can 
transit, for example, to three-mode coherence. Likewise, 
we can set up the system in some state of bimodal coher- 
ence with suitably chosen parameters K1.2 and 71,2, yet 
it can transit to higher-mode coherence. 

A formal stability criterion is only available for zero- 
mode coherence in the form of Eq. ([57]) that allows us to 
decide, without solving the EoMs, if the system is stuck 
in its alignment with B or not. 

VI. CONCLUSIONS 

We have studied two-flavor neutrino oscillations in a 
homogeneous and isotropic neutrino gas under the in- 
fiuence of neutrino-neutrino refraction. This system is 
equivalent to an ensemble of classical spins, each labeled 
by its vacuum precession frequency uj around an exter- 
nal magnetic field B and a dipole-dipole interaction with 
identical strength /i between any pair of spins. 

We have argued that the conspicuous "self-maintained 
coherence" found in this system can be understood in 
terms of an equivalent "carrier system" of a few discrete 
modes with the same dynamics. We have provided an 
explicit construction of how the original system depends 
linearly on the carrier system, assuming certain matching 
conditions. We have used this approach to construct the 
most general bimodal solution. 

Self-maintained coherence is to be contrasted with 
kinematical decoherence. We have only studied "purely 
coherent" forms of motion, but of course, depending on 
initial conditions, the system can partly or fully deco- 
here. We have not pursued the fascinating question of 
order vs. disorder in our system 6, 26]. 

Our results pertain to the simplest possible toy model 
and as such do not have any immediate practical impact 
on issues of supernova neutrino oscillations. The next 
step should be to apply these ideas to more realistic sit- 
uations, notably "multi-angle cases," where isotropy is 
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no longer assumed. The main unresolved issues in the 
theory of collective neutrino oscillations in the context of 
supernova neutrinos remains the role of multi-angle ef- 
fects. Every neutrino mode depends on energy and on 
its direction of motion, introducing much greater com- 
plications. It would be surprising if it were not possible 
to develop a more analytical understanding of collective 
flavor oscillations than has been achieved by deconstruct- 
ing numerical examples. We imagine that our results are 
only a first step toward a more complete theory of col- 
lective flavor oscillations. The ultimate goal is to provide 
a thorough theoretical underpinning for what one is ob- 
serving in numerical simulations, and perhaps eventually 
in the neutrino signal of the next galactic supernova. 



and thus the usual precession equations. 

To simplify the EoMs we introduce /3 = 5 (cji — 0^2) and 
Wc = i {uji +UJ2), leading to Pi, 2 ^ [(t^c ± (3)B + fiP] x 
Pi. 2. The only effect of ujc is a common rotation around 
B. Therefore, we go to a rotating frame where Wc = 
and thus effectively a;i_2 ~ ±13. Adding and subtracting 
these equations provides 



Q = /iPxQ and P = ;3BxQ^ 



(A4) 



where Q = Pi — P2 — {P/ IJ) B has conserved length. Up 
to a constant, the Hamiltonian becomes 



i7 = /3B-Q + ^p2. 



(A5) 
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Note Added 



As this paper went to press, a preprint appeared [28 1 
that studies in detail the properties of the quantum ver- 
sion of our Hamiltonian, Eq. and in particular its 
invariants. The constants of the motion discussed in our 
Appendix [C] are identified as the so-called Gaudin mag- 
net Hamiltonians. 



It is reminiscent of a gyroscopic pendulum with moment 
of inertia / = /i^^. The first term is the "gravitational 
potential" of the center of mass which is constrained to 
move on a sphere with radius |Q|. The second is the 
kinetic energy of the total angular momentum P. 

However, the algebraic properties of Q change. Earlier 
its Poisson brackets derived from its constituent angular 
momenta. However, the radius- vector interpretation re- 
quires [Qi,Qj\ = and [Pi,Qj\ = CijkQk, where Pi are 
the components of P. Remarkably, the EoMs of Eq. (IA4[) 
also follow with the new Poisson brackets for Q. 

In both cases, the conserved quantities are H (energy), 
B • P (angular momentum along the force direction), 
(length of the pendulum) , and Q • P (total angular mo- 
mentum projected on the radius vector or spin). The 
quantity A defined by 



(A6) 



Appendix A: Two polarization vectors 



plays the role of the natural pendulum frequency. 



We briefly review the equivalence between the dynam- 
ics of two interacting jjolarization vectors with that of a 
gyroscopic pendulum [20]. Consider the classical Hamil- 
tonian for two spin angular momenta Pi and P2 interact- 
ing with an external magnetic field and with each other 
by a dipole interaction of strength /i, 

==B-(a;iPi-|-cj2P2)+MPi-P2. (Al) 

As usual, B is a unit vector in the z-direction and a;i.2 the 
precession frequencies in the absence of /i. Each of Pi, 2 
obeys Poisson brackets which for an angular momentum 
L are [Li,Lj] = eijkLk- The Poisson brackets of the 
components of Pi with those of P2 vanish and thus P^ 2 
is conserved. Therefore, we may add {p./2) (P^ + P2) to 
the Hamiltonian and find 

if = B-(wiPi+w2P2) + |p', (A2) 

where P = Pi +P2 is the total angular momentum. The 
EoMs Pi,2 = [Pi,2,-ff] finally are 

Pl,2 = (wi,2B + ^lV)x Pi,2 (A3) 



Appendix B: Gyroscopic Pendulum 

We briefly review the textbook treatment of the sym- 
metric heavy top, also known as Lagrangian top, gyro- 
scopic pendulum, or spherical pendulum with spin. It 
is an axially symmetric body, spinning around its sym- 
metry axis (moment of inertia /a) and supported on this 
axis. Its moment of inertia relative to the point of sup- 
port is /i, mass M, gravitational acceleration g along 
the z-direction, distance I between point of support and 
center of mass, and angle relative to the z-direction. 
Therefore, V = Mg£ cos 6* is the potential energy. 

The angular momentum S along the symmetry axis 
(spin) has kinetic energy T^pin = 8^/(2/3) = I3uj'^p-^^J2. 
The point of support being on the symmetry axis pre- 
vents a torque to change 5 = | S | and so both S and Wgpin 
are conserved. The overall motion is independent of the 
internal spin angle and all tops with the same S but dif- 
ferent I3 and Wspin move in the same way. Therefore, 
we may use a single moment of inertia I = Ii = to 
describe the system. 
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The orbital angular momentum is L = / q x q, where 
q is a unit vector along the symmetry axis. It marks 
the top's orientation with zenith angle 9 and azimuthal 
angle ip. The orbital kinetic energy is 



J orbital ^2/^2 ^2 



9^ + (f^ sin^ ( 



(Bl) 



The total angular momentum J = L + S has conserved 
z-component, where Sz = S cos 9 and Lz = I(p sin^ 9. 
Here one factor of sin 9 comes from the projection of q on 
the transverse plane and the velocity is if sin 9. Therefore, 
Jz — lip sin^ 9 + S cos 9 is conserved and 



Jz — S cos 9 
I sin^ 9 



Therefore 



Torbital ^ 2 ^ ^ 



(J, -S* COS 61)2 



2/sin2 6' 

and the total energy E = T + V is 



(B2) 



(B3) 



E 



/ -o + -2 JzS cost 
2 2/sin2 6' 



Mglcos9. (B4) 



This has the form E = I9'^/2 + V{9), where V{9) is 
a potential given in terms of conserved quantities fixed 
by initial conditions. One may now solve for 9{t) by a 
quadrature and then find if{t) by integrating Eq. (jB2p . 

Next we introduce c = cos 9 as independent variable 
so that 9"^ = (? I sin^ 9 and find the usual third-order 
polynomial in c, 



o E-Mqic, 
i? ^2 (1 



Ji 



S'^-2J,Sc 



P 



(B5) 



that can be integrated using elliptic functions. A typical 
motion is nutation between two limiting angles 9i and 92 
as well as precession. When Jz ^ S = Q we have a plane 
pendulum where 



= Mg£/I 



(B6) 



transformation of an ensemble P^^ (t) obeying Eq. ([T]) . To 
motivate this transformation we observe that in the non- 
interacting case, every B • P^^ is conserved. Is there a 
similar conserved quantity in the interacting case? 

So we look for a vector X such that X-P^^ is conserved. 
A trivial example is X = Pi^, but any other X fulfilling 
X = (cjB + /iP) X X is a solution. To see this consider 
X(t) — V(t) X X(t), where V(t) is externally prescribed 
and solve for X(t) with initial condition Xq. Another 
initial condition Yq provides Y(t) and X(i) • Y{t) is 
conserved. Consider d{X-Y)/dt = X • Y -I- X • Y = 
(V X X) • Y + X • (V X Y) and after a permutation in 
one of the triple products one finds d(X • Y)/dt = 0. 

We seek X as a linear combination of the P^j and 
also of B because X = B is a solution for /i 0. 
The length of X is arbitrary, so we assume the form 
X = B + / doj' aui'^uj' and seek the coefficients a^j' ■ The 
requirement X = (wB -t- /iP) x X implies J doj' a^^'Ptj' = 
(cjB + fiP) X {B + J duj'a^'P^'). On the l.h.s. we in- 
sert Pcj' = (cj'B + /iP) X Puj', whereas on the r.h.s. 
we use B X P = B X J dut'Puj'- Collecting all terms 
we find = B X J dui' [(w — u;')atj/ — ^] Pui' and thus 
a„/ = fi/{oj - oj'). 

We are therefore led to define the linearly transformed 
polarization vectors 



P^{t) = B + /i du 



OJ ~ OJ' 



(CI) 



The integral is essentially the Hilbert transform in the uj 
variable and is understood as a Cauchy principal value. 
(The Hilbert transform of a function f{x) is defined as 
its convolution with 1/x.) So we have found a unique 
linear combination P^^ (t) of the original ensemble P^j (t) 
that obeys the original precession equation in the form 



gives us the natural frequency A. 

The motion of two interacting polarization vectors is 
equivalent to a gyroscopic pendulum (Appendix and 
therefore can be explicitly solved in terms of elliptic func- 
tions. This feat was achieved previously without invoking 
the equivalence to a pendulum p^ . 



Appendix C: Transformed ensemble of polarization 
vectors 

The functions Juj (t) derived from the carrier modes by 
Eqs. ([8|) or ([9]) represent a special case of a more general 



dtPu. = (wB + fiP) X P,^ . 



(C2) 



As a consequence, P^, P^, and P^^ • P^^ are conserved 
for any fixed fi. 

The intriguing properties of the transformed ensemble 
remain to be explored. We believe, for example, that co- 
herent oscillations correspond to Pui{t) and Pui{t) being 
coUinear at all times. 
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